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We  consider  the  problem  of  deciding  if  a  set  of  quantum  one-qudit  gates  S  =  {Ui, . . . ,  Un}  is 
universal.  We  provide  the  compact  form  criteria  leading  to  a  simple  algorithm  that  allows  deciding 
universality  of  any  given  set  of  gates  in  a  finite  number  of  steps.  Moreover,  for  a  non-universal  S 
our  criteria  indicate  what  type  of  gates  can  be  added  to  S  to  turn  it  into  a  universal  set. 


Universal  quantum  gates  play  an  important  role  in 
quantum  computing  and  quantum  optics  GS  0  [13. 
The  ability  to  effectively  manufacture  gates  operating  on 
many  modes,  using  for  example  optical  networks  that 
couple  modes  of  light  [9i,  0| ,  is  a  natural  motivation  to 
consider  the  universality  problems  not  only  for  qubits  but 
also  for  higher  dimensional  systems,  i.e.  qudits  (see  also 
[27l.  [28|  for  fermionic  linear  optics  and  quantum  metrol¬ 
ogy).  For  quantum  computing  with  qudits,  a  universal 
set  of  gates  consists  of  all  one-qudit  gates  together  with 
an  additional  two-qudit  gate  that  does  not  map  separable 
states  onto  separable  states  0  (see  UH,HB  m  for  recent 
results  in  the  context  of  universal  Hamiltonians) .  The  set 
of  all  one-qudit  gates  can  be,  however,  generated  using  a 
finite  number  of  gates  j23j-  We  say  that  one-qudit  gates 
S  =  {U\, . . . ,  Un}  C  SU (d)  are  universal  if  any  gate  from 
SU(d )  can  be  built,  with  an  arbitrary  precision,  using 
gates  from  S.  It  is  known  that  almost  all  sets  of  qudit 
gates  are  universal,  i.e.  non-universal  sets  S  of  the  given 
cardinality  are  of  measure  zero  and  can  be  characterised 
by  vanishing  of  a  finite  number  of  polynomials  in  the 
gates  entries  and  their  conjugates  [13,[2g.  Surprisingly, 
however,  these  polynomials  are  not  known  and  it  is  hard 
to  find  operationally  simple  criteria  that  decide  one-qudit 
gates  universality.  Some  special  cases  of  optical  3-mode 
gates  have  been  recently  studied  in  [5,  :32|  and  the  ap¬ 
proach  providing  an  algorithm  for  deciding  universality 
of  a  given  set  of  quantum  gates  that  can  be  implemented 
on  a  quantum  automata  has  been  proposed  [13|  (see  also 
[l ,  2,  0  for  algorithms  deciding  if  a  finitely  generated 
group  is  infinite).  The  main  obstruction  in  the  problems 
considered  in  j5|,  (33  is  the  lack  of  classification  of  finite 
subgroups  of  SU ( d )  for  d  >  4.  Nevertheless,  as  we  show 
in  this  paper  one  can  still  provide  some  reasonable  con¬ 
ditions  for  universality  of  one-qudit  gates  without  this 
knowledge. 

The  efficiency  of  universal  sets  is  typically  measured  by 
the  number  of  gates  that  are  needed  to  approximate  other 
gates  with  a  given  precision  e.  The  Solovay-Kitaev  theo¬ 
rem  states  that  all  universal  sets  are  roughly  the  same 
efficient.  More  precisely,  the  number  of  gates  needed 
to  approximate  any  gate  U  £  SU(d )  is  bounded  by 
0(logc(l/e))  0],  where  c  may  depend  only  on  d  and 
c  >  1.  Recently  there  has  been  a  bit  of  flurry  in  the 
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area  of  single  qubit  gates  |'>.  I2TI.  0, 0|  showing  that  us¬ 
ing  some  number  theoretic  results  and  conjectures  one 
can  construct  universal  sets  with  c  =  1.  The  approach 
presented  in  these  contributions  has  been  unified  in  j3l] 
where  the  author  pointed  out  the  connection  of  these  new 
results  with  the  seminal  work  about  distributing  points 
on  the  sphere  S2  [25]  that  uses  results  concerning  opti¬ 
mal  spectral  gap  for  the  averaging  operator.  Moreover, 
the  authors  of  |2Ql  showed  that  the  existence  of  the  spec¬ 
tral  gap  implies  c  =  1  for  all  symmetric  universal  sets 
of  single  qudit  gates,  where  by  symmetric  we  mean  the 
set  S  =  {U±, . . . ,  Un}  with  n  =  2k  and  Uk+i  =  U/1  for 
i  £  {1, . . . ,  k}.  Although  there  are  still  some  problems  to 
solve  in  this  area  it  seems  that  any  further  progress  would 
require  development  of  more  advanced  methods  in  pure 
mathematics  rather  than  in  quantum  information.  These 
developments  should  include  verification  of  the  spectral 
gap  conjecture.  Currently  it  is  known  to  be  true  under 
the  additional  assumption  that  gates  have  algebraic  en¬ 
tries  0|- 

In  this  paper  we  present  an  approach  that  allows  to  de¬ 
cide  universality  of  S  by  checking  the  spectra  of  the  gates 
and  solving  some  linear  equations  whose  coefficients  are 
polynomial  in  the  entries  of  the  gates  and  their  complex 
conjugates.  Moreover,  for  non-universal  S ,  our  method 
indicates  what  type  of  gates  can  be  added  to  make  S 
universal.  The  paper  is  organised  as  follows.  We  start 
from  presenting  basic  facts  concerning  the  adjoint  repre¬ 
sentation  of  SU(d).  The  adjoint  representations  assigns 
to  every  matrix  U  £  SU(d)  a  matrix  Ad;y  £  SO(d2  —  1). 
We  give  the  explicit  formula  for  Ad;/-  The  necessary 
condition  for  universality  (Lemmas  Q]  and  [5J  is  then  for¬ 
mulated  using  matrices  Ad[/  and  Ady-i,  where  U  £  S 
and  boils  down  to  checking  the  dimension  of  the  kernel 
of  the  matrix  Ms  given  by  m-  Next,  we  assume  that 
the  necessary  condition  for  universality  is  satisfied  and 
provide  sufficient  conditions  for  <  S  >  to  be  infinite  and 
thus  dense  in  SU (d).  More  precisely,  if  <  S  >  contains  at 
least  one  element  whose  Hilbert-Schmidt  distance  from 
Z(SU(d))  =  {al  :  ad  =  1}  is  both  nonzero  and  less  than 
l/y/2  then  <  S  >  is  infinite.  Combining  this  with  basic 
results  in  number  theory  we  arrive  at  our  main  results. 
In  Theorem  [T|  we  state  that  S  is  universal  if  S  contains 
at  least  one  matrix  whose  spectrum  does  not  belong  to 
some  finite  list  of  exceptional  spectra.  We  also  provide  the 
algorithm  which  allows  deciding  universality  of  any  given 
set  of  gates  S  (also  when  S  contains  matrices  with  ex¬ 
ceptional  spectra)  in  a  finite  number  of  steps.  We  discuss 
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the  correctness  of  the  algorithm  and  provide  instructive 
examples  for  S  C  SU(2). 


I.  THE  NECESSARY  CONDITION  FOR 
UNIVERSALITY 

Let  us  begin  with  introducing  the  basic  notation  used 
in  this  paper  and  explaining  the  adjoint  representation. 
The  set  of  gates  S  =  {U\, . . . ,  Un}  C  SU(d)  is  called 
universal  if  the  set  generated  by  elements  of  S 

<  S  >:=  {Ui1  •  . . .  •  Ujm  :  Uj.  £  5,  m  £  N}  (1) 

is  dense  in  SU(d),  i.e.  the  closure  <  S  >  =  SU(d).  In 
fact  <  S  >  is  always  a  Lie  group  [33].  If  this  this  group 
is  SU(d)  we  will  say  that  S  generates  SU(d). 

Let  us  denote  by  su(d)  the  Lie  algebra  of  SU ( d ).  Recall 
that  X  £  s u(d)  iff  X  is  an  antihermitian  traceless  matrix. 
Moreover,  the  Lie  algebra  s u(d)  is  a  real  vector  space 
equipped  with  a  nondegenerate  positive  inner  product 
defined  by  (X|Y)  =  —  \txXY.  For  U  £  SU(d)  and  X  £ 
su(d)  we  define 

AdtjX  :=  UXU-1. 

One  easily  checks  that  Ad;/  is  a  linear  operator  acting 
on  su (d).  It  is  also  invertible  as  (Ad;/)-1  =  Ady-i  and 
preserves  the  inner  product  as 

(Ad//A|AdyY)  =  -^trUXU-'UYU-1  = 

=  ~trXY  =  (X\Y)  (2) 

Therefore  Adjj  is  an  orthogonal  transformation  acting  on 
d2  —  1  dimensional  vector  space  su(d).  Upon  a  choice  of 
an  orthonormal  basis  {X;}“=1  in  su(d),  i.e.  basis  that 
statisfi.es  (XjjXj)  =  Sij  the  transformation  Ad;/  can  be 
expressed  in  this  basis  as  a  matrix  belonging  to  SO(d 2  — 
1),  i.e  Ad^Ady  =  I  and  det  Ad//  =  1.  The  entries  of  this 
matrix,  (Ad;/);/,  are  real  and  defined  by  the  identity: 

d 

Ad  vXj  =  U~xX3U  =  ]T(Ad//)„X;,  (3) 

i—  1 

thus  they  are  give  by 

(Adu);/  =  ~tr  (XVXp-1)  .  (4) 

Note  that  we  also  have  Ad //,//2  =  Ad;/,  Ad;/2  and  this 
way  we  obtain  the  homomorphism 

Ad  :  SU{d)  -A  SO(d2  -  1),  (5) 

that  is  known  as  the  adjoint  representation. 

For  a  set  of  d  x  d  real  matrices  M,  let  us  denote  the 
set  of  all  d  x  d  matrices  commuting  with  matrices  from 
M  by 

C(M)  =  {L  :  \L,  m]  =  0,  Vm  £  M}.  (6) 


The  adjoint  representation  of  SU(d )  is  an  absolutely  irre¬ 
ducible  real  representation  and  therefore  by  the  extended 
version  of  Schur’s  lemma  @,1H],  the  only  (d2-l)  x  (d2-l) 
matrix  that  commutes  with  all  matrices  Adgu^)  = 
{Ady  :  U  £  SU(d)}  is  proportional  to  the  identity  ma¬ 
trix,  I.  In  other  words  (^(Ads;/^))  =  {XI  :  A  £  R}. 

Example  1.  The  adjoint  representation  for  d  =  2,  i.e. 
Ad  :  SU( 2)  -a  50(3)  has  a  particularly  nice  form.  Any 
matrix  from  SU (2)  can  be  written  in  a  form 

k)  =  I  cos  (j>  +  sin  (f{kxX  +  kyY  +  kzZ),  (7) 

where  Y  =  iai,  X  =  ia2,  Z  =  ia 3  and  er,;  are  Pauli 
matrices,  k  =  ( kx,ky,kz)T  satisfies  k2  +  k2  +  k2  =  1. 
Similarly,  any  matrix  from  50(3)  has  a  form 

0((j>,  k)  =  I  +  sin  (j){—kxX12  +  kyX13  -  kzX23)+  (8) 

+2  sin2  —(—kxX\2  +  kyX\3  —  kzX23)2, 

where  Xij  =  Eij  —  Eji,  and  Eij  is  a  matrix  whose  only 
non  vanishing  entry  is  ( i,j ).  One  easily  verifies  that  the 
adjoint  representation  is  given  by 

^u(4>,k)  =  £)■  (9) 

For  U  £  SU(d),  where  d  >  2  calculation  of  matrices 
Ad;/  can  be  done  using  formula  (UJ)  upon  the  choice  of 
orthonormal  basis  in  su(d).  For  d  =  3  this  basis  is  given 
by,  for  example,  the  Gell-Mann  matrices  multiplied  by 
imaginary  unit  i.  For  higher  d  one  can  construct  an  or¬ 
thonormal  basis  of  s u(d)  in  an  analogous  way  as  for  d  =  3. 

General  considerations  that  can  be  found  in  [33j  show 
that  the  group  <  5  >  can  be  either: 

1.  <  5  >  =  SU(d),  or 

2.  <  5  >  is  infinite  and  connected,  or 

3.  <  5  >  is  infinite  and  consists  of  k  <  00  connected 
components,  where  each  component  has  the  same 
dimension  (as  a  manifold),  or 

4.  <  5  >  is  finite. 

Note  that  in  cases  1,  2,  and  3  the  group  <  5  >  has  in¬ 
finite  number  of  elements.  Thus  we  first  provide  crite¬ 
ria  that  distinguish  between  case  1  and  cases  2  and  3. 
To  this  end  we  will  use  the  adjoint  representation.  For 
S  =  {t/i, . . . ,  Un}  C  SU{d)  let  Ads  =  {Ada  :  U  £  S}. 
Note  that  if  [ L ,  Adc/J  =  0  and  [L,  Ad^2]  =  0  then 

[L,  Ad u1u2\  =  [£)  Adj/JAdpa  +  Ad;/,  [L,  Ad^2]  =  0. 

Thus  if  S  generates  SU ( d )  and  L  is  a  matrix  that  com¬ 
mutes  with  Ads  then  L  commutes  with  Ads^(d)  ■  There¬ 
fore  for  universal  S  we  have  C(Ads)  =  C^Adg/q,;))  = 
{A/  :  A  £  M}.  It  turns  out  (see  |33|)  that  the  con¬ 
verse  is  true  under  one  additional  assumption,  namely 
that  <  S  >  is  infinite. 
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Lemma  1.  For  a  set  of  special  unitary  matrices 
S  =  {Ui,...,Un}  assume  that  <  S  >  is  infinite  and 
C(Ad5)  =  {A/  :  A  £  R}.  Then  ’<S~>  =  SU(d). 

The  proof  of  this  lemma  is  based  on  the  structure  the¬ 
ory  for  semisimple  Lie  groups  and  can  be  found  in  [33j. 
Here  we  only  make  some  additional  remarks  regarding 
calculation  of  C(Ads).  Let  vec (L)  be  the  vectorisation 
of  matrix  L,  i.e.  the  vector  obtained  by  stacking  the 
columns  of  the  matrix  L  on  top  of  one  another.  One 
easily  calculates  that 

[L,  Ady]  =  0«(I®  Ady  —  Adyt  ®  I)  vec  (L)  =  0, 

where  W  is  the  complex  conjugate  and  transpose  of  U, 
i.e.  W  =  Ut.  Let 


Hilbert-Schmidt  norm  ||[/||  =  VtrUUfi  For  two  elements 
U\ ,  U2  we  have  the  following  relation  between  their  dis¬ 
tances  from  the  identity  and  the  distance  of  their  group 
commutator  from  the  identity  1 1 2|| : 

||[C/i,f/2].-/||<  V2||C/1-7||-||f/2-/||,  (13) 

[Ui,  [UltU2].],  =  I  and  ||  U2  -I  ||  <  2  =►  [Uu  U2\.  =  I. 

Let 

Ba  =  {U  £  SU(d)  :  II U  -  al\\  <  -p }  C  SU(d),  (14) 

v  2 

be  a  ball  of  radius  -fa  that  is  centred  at  elements  al.  As 

V  2 

det  al  =  1  we  need  to  assume  ad  =  1 .  Let 


Ms 


/  I  ®  Ad^  —  Ad^t  ®  I 

\  I  <8>  Adj/„  —  Ad^t  0  I 


(10) 


Lemma  2.  C(Ads)  =  {XI  :  A  £  R}  if  and  only  if  the 
kernel  of  Ms  is  one-dimensional. 


We  emphasise  the  role  of  the  adjoint  representation 
which  is  crucial  in  Lemma  [1]  In  particular  there  are  in¬ 
finite  subgroups  <  S  >  such  that  C(S)  =  C(SU(d))  but 
<  S  >  fa  SU(d).  In  Example  1  we  provide  such  a  sub¬ 
group  for  d  =  2.  We  next  characterise  space  C(Adg)  for 
SU(  2). 

Let  us  recall  that  the  composition  of  two  unitary  ma¬ 
trices  1/(7,  ki2)  =  U((j>i,  fa)U((f>2,  k2)  is  a  unitary  matrix 
with  7  and  fa2  determined  by: 

cosy  =  cos  4>  1  cos  4> 2  —  sin  fa  sin  fafa  ■  k2,  (11) 

fci2  =  — - — (fci  sin  fa  cos  (1)2  +  k2  sin  fa  coscii+  (12) 

smy 

+ki  x  k2  sin  sin  02). 

Moreover,  two  unitary  matrices  Ui(<j>i,  fci),  U2((l>2,k2) 
commute  iff  k\  ||  k2  or  f>  =  kn.  Similarly  two  orthogonal 
matrices  Oi((f>i,  ki),  02((f>2,  k2))  commute  if  k\  ||  k2  or 
one  of  (ffs  is  an  even  multiple  of  77  or  <f>i  =  ±7r  =  <f>2  and 
fci  _L  fc2.  Making  use  of  these  facts  in  j33]  we  show: 

Fact  1.  For  noncommuting  E/i(0i,fci),  U2((f>2:k2),  the 
space  C(Adt/i^i>jgi),  Adj72((^2ig2))  is  larger  than  {XI  :  A  G 

M}  if  and  only  if:  (1)  (f)  1  =  =  (j>2,  (2)  one  of  fa’s  is 

equal  to  ^  and  fa  _L  k2,  where  k  is  an  odd  integer. 


II.  WHEN  IS  <  S  >  INFINITE? 


We  next  describe  the  conditions  under  which  <  S  > 
is  infinite.  For  U\,U2  £  SU(d)  the  group  commuta¬ 
tor  is  defined  as  [Ui,U2],  =  Ufi1 .  Note  that 

[Ui,U2\  =  0  is  equivalent  to  [Ui,U2]»  =  F  The  distance 
between  elements  of  SU(d)  can  be  measured  using  the 


B=  |J  Ba.  (15) 

ad  —  l 

It  turns  out  that  noncommuting  elements  belonging  to  B 
generate  infinite  subgroups  of  SU(d ): 

Lemma  3.  Assume  that  \U\,U2}9  ^  Z(SU(d))  and 
Ui,U2  €  B .  Then  <  U\,U2  >  is  infinite. 

One  of  the  steps  in  the  proof  of  lemma  [3]  uses  relations 
m  to  show  that  the  sequence  go  =  Ui,  g±  =  [Ui,U2],, 
9k  =  [gk-i,U2],  converges  to  I  and  gn  ^  I  for  any  integer 
n  0'  (see  [33]  for  the  full  discussion). 

We  next  describe  when  U  £  Bam ,  where  am  =  e“m 
and  9m  =  =^7 r.  To  this  end  note  that 

||  U  —  amI\\2  =  2tr  I  —  otmtxU ^  —  afatrU.  (16) 

As  the  trace  of  U  is  determined  by  its  spectrum,  the 
desired  condition  can  be  expressed  in  terms  of  the  eigen¬ 
values  of  U  that  are  given  by  {e1^1, ...,  efa'1}.  fa  £  [0,  27r[ 
and  JA=i  fa  =  0  mod  27t.  Direct  calculations  lead  to: 

(17) 

i—  1 

Let  us  next  assume  that  U  £  SU (d)  does  not  belong  to 
B.  Then  one  can  show  that  there  always  exists  an  integer 
n  such  that  Un  belongs  to  some  Ba  (33j,  al  £  Z(SU(d)). 
For  a  given  U,  let  nu  be  the  smallest  integer  satisfying 
this  condition.  In  [33|  we  prove  the  modified  version  of 
the  Dirichlet’s  approximation  theorem  and  use  it  to  find 
an  upper  bound  for  Nsu{d)  :=  ma yLijnjj.  This  way,  for 
every  U  £  SU(d)  there  is  1  <  n  <  Nsu(d)  such  that 
Un  £  Ba  for  some  al  £  Z{SU{d)).  Thus  by  taking 
powers  1  <  n  <  Ngu(d)  we  can  move  every  element  of 
SU(d)  into  B.  Assume  next  that  the  necessary  condition 
for  universality  is  satisfied,  i.e.  C(Ads)  =  {XI  :  X  £  R}. 
From  Lemma  [T]  one  can  easily  deduce  that  under  the 
assumption  that  <  S  >  is  infinite  the  intersection  <  S  > 
D B  is  dense  in  B.  As  we  have  shown  in  [33}  the  necessary 
condition  for  universality  places  significant  constrains  on 
the  structure  of  <  S  >  also  in  the  case  when  <  S  >  is  a 
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finite  group,  namely  we  have  that  <  S  >  D B  is  a  subset 
of  Z(SU(d)).  Thus  <  S  >  is  finite  if  and  only  if  there 
are  no  elements  in  <  S  >  that  belong  to  B  other  than 
those  in  Z(SU(d)).  The  above  discussion  is  summarised 
by: 

Lemma  4.  Let  S  =  {Ui,...,Un}  and  assume  that 
C(Ads)  =  {XI  :  A  £  M}.  Then  <  S  >  =  SU(d)  if  and 
only  if  there  is  at  least  one  matrix  U  £<  S  >  that  belongs 
to  B\Z(SU(d)). 

We  know  that  every  element  of  S  can  be  put  to  B 
by  taking  powers  (that  are  bounded  by  Ngu(d))-  Hence 
when  <  S  >  is  finite  introducing  U  £<  S  >  to  B  must 
be  equivalent  to  introducing  it  to  Z{G).  This  condition 
can  be  phrased  in  terms  of  specra  of  the  matrices  from 
5. 


Definition  1.  Assume  U  B.  The  spectrum  of  U  is 
called  exceptional  if  it  consists  of  nth  roots  of  a  £  C 
where  ad  =  1  and  1  <  n  < 

The  set  of  exceptional  spectra  is  a  finite  set. 

To  illustrate  the  above  ideas  we  find  NSU(d\  and  the  list 
of  exceptional  spectra  for  d  =  2.  Note  that  for  any  U  £ 
SU{2)  the  spectrum  is  given  by  {e1^ ,  e~u and  therefore 
is  determined  by  one  angle  </>.  The  angle  corresponding 
to  an  exceptional  spectrum  will  be  called  an  exceptional 
angle.  Moreover,  the  centre  of  SU{2)  consists  of  two 
matrices  Z(SU{2))  =  {/,—/}.  We  start  with  recalling 
the  Dirichlet  approximation  theorem  El- 

Fact  2  (Dirichlet).  For  a  given  real  number  b  and  a  posi¬ 
tive  integer  N  there  exist  integers  1  <  n  <  N  and  p  such, 
that  nb  differs  from  p  by  at  most  >  *■  e. 

In6-pl  -  tvTT'  (1S) 


Let  [0,  2-7t)  9  (f>  =  2bn  be  the  spectral  angle  of  U .  By 
Fact  [2]  for  a  given  N  there  are  integers  p  and  1  <  n  <  N 
such  that  \nb  —  p\  <  +-j-.  Multiplying  this  inequality 

by  we  obtain  |n|  —  P%\  <  2(jv+i)  •  Note  that  (fT71) 
simplifies  to  |  sin  ^  |  <  j  or  |  sin  - 1  <  j .  Thus  for  a 
given  <f>  we  search  for  n  satisfying  |  nlf  —  Pj  |  <  arcsin^. 
Thus  <  arcsin  \  and 

2(71+1)  4 


n  < 


^  —  arcsin  | 
arcsin  \ 


=  6. 


(19) 


The  above  upper  bound  for  Nsu(d)  is  attained  for  <f>  = 
arcsin  j  (see  figure [171).  Hence  Nsu^ 2)  =  6. 

Exceptional  spectra  for  SU (2),  are  determined  by  roots 
of  1  or  —1  of  order  1  <  n  <  6,  or  equivalently,  by  primi¬ 
tive  roots  of  unity  of  order  1  <  n  <  6  and  8,  10,  12.  More 
precisely,  they  are  given  by  {el9 ,  e~10},  where 


9  £  {kn, 


/C27T  k^TT  kd7T  ksTT  k^K , 


(20) 


FIG.  1.  Condition  (1171)  for  SU( 2).  The  gray  segments 
are  determined  by  |  sin  f|  <  I  and  the  white  segments  by 

I  sin  VI  <i 


and  gcd(fc+)  =  1.  The  number  of  exceptional  angles 
can  be  calculated  usinjr  the  Euler  totient  function,  and 
is  equal  £T=1  <t>{i)  +  X]j=4  <^(2*)  =  24.  For  higher  dimen¬ 
sional  groups,  as  we  discuss  in  (33|,  the  number  Ngjj ^ 
grows  exponentially  with  d.  Our  main  result  is: 

Theorem  1.  Assume  that  C  (Ad^ , Adc/n)  =  {A  I  : 
A  £  R}  and  at  least  one  matrix  Ui  has  a  nonexceptional 
spectrum.  Then  <  Hi, . . . ,  Un  >  =  SU (d). 

III.  THE  ALGORITHM  FOR  DECIDING 
UNIVERSALITY  OF  S  =  {Ui,  ...,[/„}  C  SU(d) 

The  case  when  all  matrices  {Ui, . . . ,  Un}  have  ex¬ 
ceptional  spectra  requires  an  algorithm  which  we  next 
present.  Our  algorithm  allows  deciding  universality  of 
any  given  set  of  SU(d)  gates  in  a  finite  number  of  steps. 

The  algorithm 

Step  1:  Check  if  C(Ad^)  =  {XI  :  A  £  R}.  This  can 
be  done  by  checking  the  dimension  of  the  kernel  of 
the  matrix  Ms  cud  constructed  from  the  entries 
of  matrices  {Adj/i; . . . ,  Ad[/n}  and  thus  is  a  linear 
algebra  problem.  If  the  answer  is  NO  stop  as  the 
set  S  is  not  universal.  If  YES,  set  l  =  1  and  go  to 
step  2. 

Step  2  :  Check  if  there  is  a  matrix  U  £  S  for  which  Unu 
belongs  to  B  but  not  to  Z(SU(d)),  where  1  <  nu  < 
Nsu(d)  ■  This  can  be  done  using  formula  (fl6l) .  If  the 
answer  is  YES  S  is  universal.  If  the  answer  is  NO, 
set  l  =  l  +  1. 

Step  3  :  Define  the  new  set  S  by  adding  to  S  words  of 
length  Z,  i.e  products  of  elements  from  S  of  length 
l.  If  the  new  S  is  equal  to  the  old  one,  the  group 
<  S  >  is  finite.  Otherwise  go  to  step  2. 

The  major  advantage  of  our  approach  is  the  fact  that 
we  can  make  decisions  in  steps  2  and  3  in  finite  ‘time’. 
It  is  also  clear  that  for  randomly  chosen  matrices  S  = 
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{U\, . . . ,  Un}  C  SU(d)  our  algorithm  terminates  with 
probability  1  in  Step  2  for  1  =  1.  This  is  a  direct  conse¬ 
quence  of  the  fact  that  exceptional  spectra  form  a  finite 
set. 

Let  us  next  discuss  the  correctness  of  our  algorithm. 
Assume  that  S  passes  positively  the  necessary  condition 
for  universality,  i.e.  the  Step  1.  If  the  group  <  <S  >  is 
finite  the  algorithm  terminates  in  Step  3  for  some  finite 
l.  On  the  other  hand,  as  a  direct  consequence  of  Lemma 
SI  for  an  infinite  <  S  ^  the  algorithm  must  terminate  m 
Step  2  for  some  finite  l.  One  can  also  argue  that  if  all 
finite-length  words  have  exceptional  spectra  then  they 
cannot  form  a  dense  subset.  Thus  if  <  S  >  is  dense 
then  it  must  contain  words  of  a  finite  length  that  have 
non-exceptional  spectra.  These  words  terminate  the  al¬ 
gorithm  in  Step  2.  Moreover,  we  have  the  following: 

Fact  3.  Assume  <  S  >  is  dense  in  SU(d).  The  length 
of  a  word  that  gives  termination  of  the  universality  al¬ 
gorithm  is  at  most  the  length  l  such  that  words  of  length 
k  <  l  form  an  e-net  that  covers  SU(d),  where  e  =  2v/g+g 
and  6  >  0  is  arbitrary  small. 


both  <j)i  and  02  are  exceptional  angles  and  the  product 
U((f>i2,  kif)  =  t/(0i,fci)t/(02,fc2)  lias  exceptional  0  12  as 
otherwise  the  algorithm  terminates  in  Step  2  with  either 
l  =  1  or  l  =  2.  For  each  such  pair  (there  are  finitely  many 
of  them)  we  find  the  Step  and  l  that  gives  termination 
of  our  algorithm.  The  detailed  discussion  of  the  results 
and  their  connection  to  finite  subgroups  of  SU( 2)  can 
be  found  in  [33].  In  short  words,  it  turns  out  that  our 
algorithm  terminates  at  Step  2  if  and  only  if  1  <  l  <  4 
and  at  Step  3  when  5  <  l  <  13. 

Fact  4.  For  S  =  { t/ (<^>i ,  fci ) ,  C7 (02 ,  ^2)}  C  SU{2)  the 
algorithm  for  checking  universality  terminates  for  l  < 
13.  Moreover,  the  set  S  is  universal  if  and  only  if  the 
algorithm  terminates  for  l  <  4. 

The  main  conclusion  from  Fact  H]  is  that  one  can  de¬ 
cide  universality  of  any  two-element  subset  of  SU (2)  by 
looking  at  words  of  the  length  at  most  4. 

IV.  EXAMPLES  FOR  SU( 2) 


Proof.  Assume  that  words  of  the  length  k  <  l  built  from 
elements  S  form  an  e-net  for  SU(d),  where  e  =  2v/g+a 
and  S  >  0  is  arbitrary  small.  Let  U  be  an  element  of 
SU(d)  whose  distance  from  the  identity  is  exactly 
(see  Figure[5]).  Then  by  the  definition  of  e-net  there  must 
be  at  least  one  word  w  £<  S  >  of  length  k  <  l  contained 


in  the  ball  C  of  radius  e  = 


centred  at  U.  But  this 


2V2+S 

ball  is  contained  in  B\  \  I.  Hence  w  gives  termination  of 
the  universality  algorithm  in  Step  2.  The  result  follows. 

□ 


In  the  remaining  part  of  this  paper  we  demonstrate  our 
approach  calculating  a  few  examples.  They  are  chosen 
particularly  to  elucidate  the  importance  of  the  conditions 
given  by  Theorem  [Tj 

a.  Example  1  Let  S  =  {(7(0,  k\),  U(n/2,  /C2)}, 
where  k\  _L  and  <f>  1  is  an  irrational  multiple  7 r.  For 
example,  when  k\  =  (0,  0, 1)  and  k2  =  (1, 0,  0),  we  have 

U(cf,  h)  =  (  ^  )  ,  U{ it/2,  k2)  =  (  _°j  J  )  - 

(21) 


FIG.  2.  The  proof  of  Fact  [3] 


The  formulation  of  Fact  [3]  is  related  to  the  results  con¬ 
tained  in  [l8|. 

In  order  to  demonstrate  how  efficient  is  our  algorithm 
we  determine  the  maximal  l  which  gives  its  termination  in 
Step  2  and  Step  3  respectively  for  SU(2).  For  simplicity 
we  consider  S  of  the  form  S  =  {U((f>i,  ki),  U{ej>2,  £2)}  C 
SU (2).  To  this  end  it  is  enough  to  consider  the  case  when 


U((j>,  ki)  is  of  an  infinite  order  and  since  U((j>,  k\)  and 
U (7t/2,  /C2)  do  not  commute  we  have  that  <  S  >  is  infinite 
and  not  abelian.  By  Fact  [T|  however, 

^(Adj^  ^),  Ady^/^jk))  7^  '■  ^  S  C}, 

and  hence  <  S  >  SU( 2).  For  example  0(7r,  ki)  £ 
SO(  3)  commutes  with  both  Ad[/^i  £,  =  0(201,  hi)  and 
Adf7(7r/2,fc2)  =  0(TT,k2).  Interestingly,  however 

C(U1,U2)  =  {XI:\eC}. 

To  understand  the  structure  of  the  group  <  S  >  note  that 
U(n/2,  ^2)17(01,  ki)U~1(n/2,  fa)  =  Hence 

U(n/2,k2)  is  a  normaliser  of  <  t/(0i,/ci)  >.  Thus 
the  group  <  S  >  consists  of  two  connected  components. 
The  first  one  is  given  by  one- parameter  group  U[t,k\ ), 
where  f  €  1  and  the  other  one  by  elements  of  the  form 
f7(7r/2,  fc2)I7(t,  ki).  The  adjoint  representation  is  able 
to  identify  infinite  disconnected  subgroups  whereas  the 
defining  representation  is  not.  Moreover,  we  know  ex¬ 
actly  how  to  fix  non-universality  of  the  set  S.  For  exam¬ 
ple,  we  can  add  one  matrix  1/(7,  fc7)  such  that  7  kir 
and  fc7  is  neither  parallel  nor  orthogonal  to  k\  and  &2 . 
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b.  Example  2  Let  H  be  the  Hadamard  gate  and 
X^  a  phase  gate  with  an  arbitrary  phase  <f>: 


1  1 
1  -1 


e-^  0  \ 

0  J  ' 


(22) 


Using  our  notation  H  = 

^(0,1,1)  and  X0  =  U{</>, 

ki  ■  = 


i 

YU 


U(n/2,kH ),  where  fc#  = 
fcTJ,  where  kT</i  =  (0,0,1), 


Our  goal  is  to  check  for  which  </>,  <  Id,  X^  >  =  SC/ (2). 


non-exceptional  spectrum  and  whose  fc  is  neither 
parallel  nor  orthogonal  to  kn  and  kpn/2. 

Case  3:  For  <f>  ^  again  by  Fact[0  C(Ad/r,  AdT(^))  = 
{XI  :  A  £  R}  and  we  just  need  to  check  if  <  H,  > 
is  infinite.  We  distinguish  three  possibilities: 

1 .  We  first  assume  that  <f>  is  not  exceptional.  Then 
by  Theorem  Q]  (H,  X^)  =  SU(2).  Our  algorithm 
for  deciding  universality  terminates  at  step  2  with 
l  =  1. 

2.  We  next  consider  the  exceptional  angles.  For 


ffc3  7T  fc57T  k6TT 
e{~’~  ~}’ 


gcd  (fcj,i)  =  1, 


we  look  at  the  product  U('y,kHr)  =  HT^  = 
U(tt/2,  ku)U((j>,  fcr)-  Using  formula  (fill)  we  calcu¬ 
late  cosy,  compare  it  with  cos ip  for  all  exceptional 
angles  ip  and  find  out  they  never  agree.  Hence 
7  is  not  exceptional.  Thus  by  Theorem  [T|  we  get 
<  HTfj,  >  =  SU{ 2).  Our  algorithm  for  deciding 
universality  terminates  in  Step  2  with  1  =  2. 

3.  We  are  left  with  <f>  =  where  gcd(&4,4)  =  1. 
There  are  exactly  four  such  angles.  Calculations  of 
U (7,  k hTj,  )  =  HT^  shows  that  7  is  exceptional,  i.e 
7  =  where  gcd(fc3,3)  =  1.  Moreover,  taking 
further  products  results  in  a  finite  subgroup  con¬ 
sisting  of  48  elements  (all  have  exceptional  spectra) 
known  as  the  binary  octahedral  group.  Our  algo¬ 
rithm  for  deciding  universality  terminates  in  Step  3 
with  l  =  8.  Fixing  non-universality  can  be  accom¬ 
plished  by,  for  example,  adding  one  gate  U{ip,k^) 
with  a  non-exceptional  ip  and  an  arbitrary  k 


Case  1:  If  <p  =  kir  then  X^  =  ±1  and  the  generated 
group  is  the  finite  cyclic  group  of  the  order  4  when 
(f>  =  0  or  the  order  8  when  <p  =  7 r. 


As  we  can  see  from  the  above  example  our  algorithm 
requires  at  most  words  of  length  l  =  8  to  terminate  for 
any  H  and  X^. 


Case  2:  When  cp  =  ^  and  k  is  odd,  by  Fact|T|  we  have 
that  C(Ad#,Ad7y)  is  larger  than  {A/  :  A  £  R} 
and  hence  <  id,  Xs^ L  >  ^  SU(2).  In  fact  it  is  the 
finite  dicyclic  group  of  order  16  whose  generators 
are  HT  and  X.  Fixing  universality  in  this  case 
requires,  for  example,  adding  a  matrix  that  has  a 
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